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Abstract 



We present a computation of the k-string tension m the large N limit of the two- 
dimensional lattice Yang-Mills theory. It is well known that the problems of computing 
the partition function and the Wilson loop can be both reduced to a unitary matrix 
integral which has a third order phase transition separating weak and strong coupling. 
We give an explicit computation of the interaction energy for fc-strings in the large 
N limit when -fe is held constant and non-zero. In this limit, the interaction energy 
is finite and attractive. We show that, in the strong coupling phase, the k — )■ N — k 
duality is realized as a first order phase transition. We also show that the lattice 
/c-string tension reduces to the expected Casimir scaling in the continuum limit. 

(N 
> 

■ K-string tension has been proposed as an interesting probe of the confining phase of Yang- 
Mills theory [1] [2] [3] [3] [5] [6] [T] [8j [9j [1Q] . In the dual superconductor picture of confinement, 

■ lines of electric flux which emanate from a source with color charge are confined to flux 
tubes, the confining strings. If a colored source has center-charge k, it is the endpoint of 
k flux tubes. There is immediately an interesting question as to whether these flux tubes 
attract and perhaps bind together to form a single tube with k units of flux, or whether 
they repel each other and tend to remain as individual vortices, that is whether the dual 

^ ■ superconductor is of type I or type II, respectively. 

In Yang-Mills theory, or any gauge theory with only adjoint (or other center-neutral) 
matter, the k-string tension au is defined as 

e -c k A = (Tr Ak e^o A ^ 
dimA k \ h 

where C is typically a rectangular loop with dimensions LxT, A = LT is the area subtended 
by that loop and Ak is the irreducible representation of the gauge group where the Young 
tableau consists of a single column of k boxes, depicted in Fig. [TJ We are assuming that the 
gauge group is SU (N). The physical interpretation of <Q lies in the idea that, if we introduce 
a heavy fc-quark and a heavy anti-/c-quark separated a distance L much greater than the 
confinement scale into the confining phase of Yang-Mills theory, the energy attributable to 
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Figure 1: The Young tableau corresponding to a completely antisymmetric representation 
Ak of SU(N) consists of a single column with k boxes. 

the gauge theory mediated interaction between them is E = <JkL. The linear growth of 
the energy with distance is a signal of confinement. The coefficient at is a measure of the 
strength of the confining interaction. 

An intuitive way of thinking about the fc-string is to imagine that we bring k confining 
strings close to each other. If these elementary strings have an attractive interaction, er& < 
ka±, and the strings can from a bound state which we call the fc-string. The fc-string tension 
is not easily accessible to perturbation theory and, aside from some supersymmetric models 
where exact results can be obtained [TT], it is most readily studied by numerical simulations 
of lattice gauge theory. In this Letter, we will study it in a solvable model, the large N 
limit of two-dimensional Yang-Mills theory that was originally formulated and solved by 
Gross and Witten [12J. The fundamental and adjoint Wilson loops in this model have 
been used as models of deconfinement transitions in finite temperature Yang-Mills theory 
Khokhlachev: 1979tx[H] [B] . 

If we were to consider an expectation value of a Wilson loop as in ([Q), but taken in 
some other representation of the gauge group where the Young tableau also has k boxes, 
the expectation value should turn out to be independent of the representation and equal 
to that obtained for a k-string in the completely antisymmetric representation Ak- This is 
because we expect the color charge of the quarks to be screened by the gluons which, since 
they are center charge neutral, conserve the center charge k but, by combining with the 
strings of electric flux, can alter the precise representation. This screening is not expected 
to occur in two dimensions where there are no dynamical gluons. It is also not expected to 
occur in the large N limit where large N factorization means that mixing of representations 
is suppressed by factors of j^. Lack of mixing of representations would lead us to expect 
that in two dimensions, or at large N, the k-string tension is always given by the weak 
coupling result, ak = k<j\. At a first glance, large N factorization would suggest that the 
interaction energy of a fc-string, ak — ka±, should be of order in the 't Hooft large N 
expansion, and therefore vanishingly small in the large N limit. However, since at infinite 
N all representations have the same tension, lifting the degeneracy at finite N can lead 
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to corrections of order 1/N, as we will see shortly 0. There are lattice simulations which 
suggest that corrections to the fc-string tension in 3 dimensions indeed scale like rather 
than 4j [TD]. In addition, the lattice simulations seem to indicate that the splitting of the 
energies of representations with the same center charge is identifiable, is also of order -k and 
is unexpectedly small. 

The /c-string tension is expected to be symmetric under the replacement of k by iV — k so 
that 

0"fc = &N-k ■ (2) 

This is obtained by replacing probe quarks by anti-quarks and vice- versa. In contrast, weak 
coupling implies that a k ~ ka%. While this fundamental k — >■ iV — k symmetry is compatible 
with /c-strings being weakly coupled for k ^ N, since k and N — k cannot both be much less 
than N, in the regime where N is large and k ~ N, /c-strings cannot be weakly interacting. 

Continuum Yang-Mills theory in two dimensions is solvable and the fc-string tension is 
known exactly, 

of = 2 = \C 2 (U k ) , (3) 

where A = g 2 N is the 't Hooft coupling and g is the Yang-Mills theory coupling constant. 
C^-R) is the quadratic Casimir operator of the representation TZ k - As expected, because 
there are no dynamical gluons, the string tension is not independent of representation. For 
a completely antisymmetric representation A k , it is 

CMk) = ^ • (4) 

This expression is symmetric under the replacement k — >■ N — k. This fc-string tension has 
a® =2 = ka® =2 — ■^■crf )=2 , the first term interpreted as the energy of k non-interacting fun- 
damental strings and the second term being the interaction energy. Thus, the interaction 



1 In the infinite N 't Hooft limit, and with small values of k, all representations with the same center 
charge should have the same string tension, ak = ka±, that is, they are degenerate. These representations 
are mixed with each other and the degeneracy is resolved by non-planar corrections where the "interaction" 
is of order . A simple model for this mixing is the two-state system where the Hamiltonian is 



°) + f0 9/N 



E 2 J \g/N 



Corrections to the energy are give 



Ei — > E\ — 777-^ — I" ■ ■ • 1 E% — > E2 



N 2 E-2 - Ei " ' ' N 2 E 2 ~Ei 

However, if the leading order is degenerate, E\ = E% = Eq then 

Ei -> E - — + . . . , E 2 ->• E Q + — + . . . 

Generically, corrections are of order 1/N. (A version of this argument appears in Ref. [8].) We thank Barak 
Bringoltz for comments which motivated us to clarify this argument. 
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between /c-strings in two dimensions seems to be attractive when they are in the antisym- 
metric representation. It is also easy to see that for expectation values of a Wilson loop in 
other representations, the interaction energy can be repulsive. 

The Casimir scaling hypothesis is that the fc-string tension in higher dimensions remains 
proportional to C^GAjfe), dependent on the center charge k but, because of screening, oth- 
erwise independent of the representation of the probe quarks. This hypothesis is seemingly 
contradicted by the sine-law: an alternative, first principles result found in softly broken 
M = 2 supersymmetric Yang-Mills theory [IT] where the fc-string is a BPS object 



0~k 



mAN sin . (5) 



The sine-law is consistent with results found in MQCD [16] and AdS/CFT p2j. One inter- 
esting difference between the Casimir scaling and sine-law behaviors is that the corrections 
to the leading order in the large N limit of Casimir scaling are of order 1/N whereas for the 
sine-law they are of order 1/N 2 . In both cases, though, k-strings are strongly interacting 
when k ~ iV as was discussed above. 

In the following, we shall find a third expression for the k string tension which follows 
from the large N limit of two dimensional lattice Yang-Mills theory. Above, we have pointed 
out two non-generic features of the fc-string in two dimensions and at large N. One is the 
absence of dynamical gluons in two dimensions which leads to representation dependence, 
as is already seen, for example, in the continuum formula 0- The other is the absence of 
mixing of representations which happens in the large N limit in any dimensions, so that 
/c-strings are non-interacting. In the following, we shall find a way to circumvent the second 
shortcoming by studying representations with center charge k where k ~ N. As has been 
discussed, this is not compatible with weak coupling; mixing of representations, which would 
be suppressed by powers of 4, is allowed and is of unit magnitude when k ~ iV '. Since the 
interaction is not suppressed by a power of jj, one question that our computation can answer 
is whether confining strings continue to have an attractive interaction so that fc-strings are 
stable when the interaction is strong. We shall find that it is indeed attractive and — ko\ 
is negative for all values < < 1. We shall also confirm that the string tension in the two- 
dimensional model indeed depends on the representation and, for the completely symmetric 
representation with center charge k, the interaction between confining strings is repulsive for 
all values of j? > 0. 

A Wilson loop in the two dimensional lattice Yang-Mills theory is computed by the integral 
Wn[c] = I Ef Tr(n ^ Ut+U ^ P Ul) Ueec Uj (g) 

where P and t label plaquettes and links, respectively, of a square lattice, Ue is a unitary 
matrix residing on the links, [dU(\ is the invariant Haar integration measure for SU(N), TZ 
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is an irreducible representation and Tr^ is the character in that representation. This model 
was first solved in the large N limit by Gross and Witten [T2]. Using the gauge fixing of 
Ref. [123 and the Peter- Weyl theorem J[dU]U%Ulf = , it is possible to reduce this 

to a one-plaquette unitary matrix model [18] 



(7) 



W n [C] = [w n [C]} A = e 

N 



o n = - In 



dim7£ 

j[dU}e^< u+W ) 



where A is the area of the loop - the number of plaquettes in an area that is bounded by 
the loop - and the string tension a-ji is measured in units of inverse lattice spacing squared. 
Note that the Wilson loop in the one plaquette model is normalized by the dimension of the 
representation. This is a result of repeated use of the Peter- Weyl theorem when reducing to 
the one-plaquette model. It has the important consequence that, since in this simple one- 
matrix model, the dimension of the representation is the upper bound on the expectation 
value of the character of a matrix in that representation, the string tension is always positive 
and the gauge theory interaction is always confining. 

The one-plaquette model which computes a-ji in (j7j) is solvable in the large N limit. 
Gauge invariance allows one to diagonalize U = diag (e 1 ^ 1 , . . . , e^ N j and the remaining 
integral over eigenvalues can be done in the saddle-point approximation where the resulting 
action is of order N 2 , so that fluctuations are suppressed by powers of j^. The classical 
configuration of the eigenvalues is determined by minimizing the action plus a term arising 
from the integration measure. In all of the cases that we shall consider, if we seek only the 
leading order in large N, the loop inserted into (j7|) can be treated as a probe where the 
expectation value is evaluated simply by plugging in the classical distribution of eigenvalues 
that is determined by the action. This distribution is characterized by the eigenvalue density, 
defined as the large iV limit of p((p) = i J2i &{4> ~ 4>i)- The model (J7|) has a third order phase 
transition when A = 2. The two phases have eigenvalue densities 

p{4>) = ^(l + f cos0) , for A>2 , 



p{4>) = ^cosf^-sin^f , forO<A<2, (9) 

where, in (jHJ), |</>| < 2 arcsin and on the remainder of the circle p = 0. 

The easiest representations to consider are the totally antisymmetric Ak whose Young 
tableau consists of a single column of k boxes and the totally symmetric Sk whose Young 
tableau is a row with k boxes. The analysis of the symmetric representation in a unitary 
one-matrix model was explained in detail in Ref. [19J and we will refer the reader there for 
the details. Here, we will concentrate on the antisymmetric representation, corresponding 
to a k-string. 
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In terms of eigenvalues, the symmetry can be used to order the indices in the trace so 
that they are non-decreasing 



e t<pa 1 e t<pa 2 e l( Po. k 



ai<—<a k 

It is convenient to obtain this expressions from a generating function [20j[19j 

N 



(10) 



J a=l 



where the contour integral projects onto the term in a Taylor expansion of the integrand 
which contains k eigenvalues. The covariant expression is 



(TrA k U) = y^ T (exp[Trln(l + t[/)]). 
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In the large N limit, the traces in the exponents in the above equations are replaced by 
integrals over the eigenvalue densities (JSJ) or ([? 



(^A k U) 



dt 



exp 



N / d(pp((p) ln(l + te**) — hint 



(16) 



2 The reader might have the concern that the presence of the loop variable in the path integral, though 
it does not alter the eigenvalue distribution to the leading order TV , it will have an effect at order 1/N and 
a 1/JV correction in the order N 2 part of the action would contribute a term of order N which competes 
with the string tension which we are computing. To see why this is not a problem, consider the tension in 
the large N limit is given by 



a k = inf 

(p,t) 



N 2 S[p] + A J p-X-N J pin (1 + te^) + hint 



-inf 

p 



N 2 S[p] + A / p - A 



(13) 



where S[p] is the effective action consisting of S s plus a contribution from the integration measure. The 
saddle-point equations are 



SS_]_ 
Ip~~N 

for the first infimum and 



In (1 + te*+) + 



N 2 







1 



' 1 + te** 



k 1 



5S 
5p 



A 

N 2 



(14) 



(15) 



for the second infimum. The eigenvalue density which satisfies (II 5[) is po. Then the density which satisfies 
([T4"l) differs from it by a correction of order -A., po + j?Pi- However, since po satisfies (fT5"]) . it is easy to see 
that if we are interested in <t/- only to accuracy of order N, we can simply use po m the equation which 
determines t and, to the same accuracy (where we trust the order N but not the order N° contribution), in 
the expression for cr^ in (|13|) . This justifies our use of the "probe approximation" where we use the eigenvalue 
distribution of the effective unitary matrix model to compute the generating function. We note that a similar 
probe approximation is made when analyzing the dual objects on the string theory side of the AdS/CFT 
correspondence in Ref.|19j. 
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When k is large, we can use the saddle-point approximation to evaluate the integral over t 
in (I12p . Let t satisfy the saddle-point equation 



d(f> P U) ^-r = T7 • (17) 

The functions R^ k (t) in (I17p are related to the resolvent of the matrix model and are holo- 
morphic functions of t with cut singularities on the unit circle determined by the support of 
p(4>). With the strong coupling phase (jSJ), A > 2, the cut singularity on the left-hand-side of 
(fT7|) occupies the entire unit circle and divides the plane into two regions, the interior and the 
exterior of the unit circle. We extend each of these regions to disjoint Riemann sheets, each 
of which covers the entire plane. We then search for solutions of the saddle-point equations 
on each of these sheets. The saddle-point equations for each of the two sheets are 

¥ = *> (18) 



il k_ 

A 2 N 



(19) 



and the solutions are 



h = Xjf , (20) 



N ' 

*2 = \l^kjN ' ( 21 ) 

respectively. Using these we find two branches of the string tension, each coming from 
anti-derivatives by t of (JI8J) and f|T9|) respectively. The possibilities are 

o k = -k - (N - k) ln(l - ±) + k In A and a k = - (N - k) + k In £ + (iV - fc) In A . (22) 

In principle, to get the correct result, we must choose the one that has the smaller value of 
the tension, opening the possibility of a phase transition. However, we need to be careful 
because using a smooth density function p in f fl6|) might lead to the wrong result when k/N is 
large enough. A smooth p(<p) is an approximation to A^ 1 5(<f) — <p a ). To accurately obtain 
a k th Taylor expansion coefficient of t in the argument of the exponent in equation (|16|) . we 
need to know the k th Fourier coefficient of <fi in p{4>). When k/N (which we are asuming to 
be of order 1 as A^ — > oo) is larger than the inverse of the density of eigenvalues in some 
region, replacing the discrete sum over eigenvalues with the continuum density might not 
lead to the correct answer, since the wavelength of the fourier mode we are interested in will 
be shorter than the inter-eigenvalue spacing. Equation (fT6]) thus should not be trusted for 
k/N large enough. However, if we analiticaly continue the t-integral onto the second sheet, 
we are effectively exchanging k and N — k, and the second sheet answers are trustworthy for 
small (N - k) /N. 
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To test whether the results in equation ([221) are in agreement with a discrete eigenvalue 
density, we have evaluated the expression in equation (ITSj) for = N^ 1 Xla'K^ - 0a) "with 
<p a taken to match the large- N distribution (jSJ) at some finite N. Fig. |2] shows the results. 
For larger values of A (A > e = 2.718 . . .), the first sheet result dominates for k/N < 1/2 and 
the second sheet result dominates for k/N > 1/2, with a phase transition at k/N = 1/2. 
The existence of the phase transition is well supported by the computation with a discrete 
density, especially at larger values at A where the phase transition is sharper. For 2 < A < e, 
it would appear that there are two more phase transitions as the saddle points exchange 
positions again. It is clear from the discrete results, however, that these are an artifact 
of the continuum approximation. As described above, we cannot trust the green curve for 
larger values of k/N, or the blue curve for smaller values of k/N. 

In summary, cr^ depends on whether is greater than or less than ~, 



These are accurate to the leading order in large N, 0(N), when k ~ N and we expect 
that there are corrections of order one, O(N ). The result has an explicit k — >• N — k duality 
which is obtained by exchanging the saddle points. This exchange of saddle points is a first 
order phase transition (the first derivative by is discontinuous at 4 = |). 

We recall that for a single box fundamental Wilson loop in the strong coupling phase of 
the Gross- Witten model [12] . G\ = In A. In the top expression in (1231) . we can clearly identify 
the contribution from the energies of the k fundamental strings, ko\ = A; In A. We can also 
confirm from the remaining contributions in ( )23|) that the /c-string tension is indeed always 
less, that is < kai, as expected for a stable fc-string. In fact, for small jj, 



As in the continuum string tension, the correction is of order 1/N, rather than 1/N 2 . 

It is interesting to compare the string tensions in ( |23l) with those in the symmetric repre- 
sentation [19] 



The symmetric string tension is always larger than the tension of k single strings, as k > 
ka\ > cTfc, even for small 4 where as k = ko\ + |4 + . . .. This is consistent with the idea 
that representations other than the antisymmetric one should have higher energies. Also, as 
expected, the symmetric representation tension is not identical to the antisymmetric one, as 
there are no dynamical gluons to implement screening. In this representation, the strings 
seem to repel each other. 




(23) 



(24) 



a Sk = (N + k)\n(l + ±) + kln\-k , A<2 . 



(25) 
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X = 8 



X = 4 




X = 2.5 1 = 2.2 




Figure 2: as a function of k for iV = 500. The green and blue solid lines are obtained 
from the first and second terms in equation ([221 respectively, while discrete points (black) 
are computed directly from equation ( JIB]) with a discrete density p(4>). 
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Now let us consider the weakly coupled, gapped phase with eigenvalue distribution (J5J). 
The string tension is obtained by integrating ( II 7p to get the exponent 



d<f>p{<i>) In (1 + ie**) + £ln* = (26) 
(* + !) (i+l-yf(i+l)*-2\t) i+l-sj{i+lf-2\i ! k 



In - 1 h — In t 



2Xt Xt 2 N 

where t is adjusted to be an extremum of this expression. Setting the derivative of the above 
expression by t to zero yields an equation that is identical to the result of taking the integral 
over the eigenvalue density (jHD in the saddle-point equation ( flTj) . 



i + 1 - yj(i+ I) 2 -2\i 
~2~\t 



#-! = (*-!) ^"7777 ' , A<2. (27) 



The position of the cut singularity in the square root is determined by the cut singularity of 
the resolvent in the saddle-point equation (fTTl) . The cut is on the unit circle, antipodal to 
the region where the eigenvalue density has support. The cut crosses the negative real axis 
at t — —1. The sign of the square root has been adjusted so that it has the right behavior at 
t ~ 0. As we follow the negative real axis, the sign must flip at t = — 1 so that the expression 
has the correct large (negative) t behavior. There are two solutions 

1 - - (l - 2—) 2 1 - 2— 
1 4 V 1 z n) _i_ n 



4 = 4 / ± ; " V£ (1 - 2^Y + 1 - | . (28) 

2(1- |) 2(l-^) Vl6V N) 2 V ' 

The argument of the square root in this expression is positive in the domain of parameters 
of interest. Both roots are non- negative real numbers which can vanish only when j? = or 



N 

jj = 1. They are related by 

M- = T^T >*-(*) = I fc \ • (29) 

7v l + V ~ n) 

Note that this is similar to the relationship in (120]) and (12 ip. For real values of the parameters 
A and the argument of the square root in (j2SJ) is always positive, so this square root is 
unambiguously defined in the entire region of interest. It does have two branches which can 
both be confirmed to satisfy the saddle-point equation. It's the smaller root, t_, which lives 
on the first sheet and goes to zero as 4 — > 0, that will turn out to be the correct one to use 
based on the discrete density calculation. 

We can substitute £_ into (1281) . taking into account the normalization by adding In diuuA/c = 
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Figure 3: as a function of k for N = 500 in the weakly coupled phase. The green and blue 
solid lines correspond to equation [30] and a similar expression obtained for the other sheet 
of the saddle point equation. The green line represents our result for the /c-string tension 
in the semi-circle distribution. Discrete points (black) are computed directly from equation 
( TT6l) using a discrete eigenvalue density. 



ln(N\/k\(N — k)\) to the result, to get the tension 



N 



-Ai n A_(i_A) ln (i_A)_ 



N 



16 



(1 



Z N ) 



+ In 



A 



l + f(l-2 fc 



+ 



k 1 

iv ln - 



N - 



Z N ) 



2 ~ 16 ^ 



oA^ 



N> 



2 ~ 16 ' 



1-2A) 2 



2(1 



A" 

7V' 



(30) 



Similar expression can be obtained starting with the other root, t + . Fig. ([3]) summarizes 
the results: even though the second sheet saddle point is smaller, it is the first sheet result 
stated above which agrees with a discrete eigenvalue density computation. 

Though it is not immediately obvious from expression ( 1301) , the tension is symmetric under 



replacement of 4 by 1 



•4. This is clearly seen in Fig. [3j 



For small 



<jfc = k In 



1 



fc 2 A(8-3A) 



1-1 



(31) 



iV 2(4 -A) 2 

where, as expected, the string tension for a k — 1 string in the weak coupling phase is 
<7i = In 



i 



— j- [12] . As in the strong coupling phase, the interaction term is attractive and is 
1_ I 
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suppressed by one factor of 1/N. We can also confirm that the interaction is attractive for 
all values of < < 1 and for all < A < 2. 

On the other hand, in the weak coupling phase, the symmetric representation string 
tension can be deduced from the results in Ref. [T9l 



1 



+ 



.il n i + 

N N 



;i+ *)in(i + £) + - 



i + t(i 



i-t- 



2^) 



1-f + fi(l + 2A)2 



N 



In 



-l + 4(l + 2^) 2 + (l 



N ' 



2(1 + §) 



' N' 



2#)\/l-| + fi(l + 2#) 2 



(32) 



By considering the small limit, <xs fc 



fcln 



+ 



fc2 A(8-3A) 



+ o 



we see that, in 



N 2(4-A) T 

the symmetric representation, the interaction is repulsive. It is possible to show that it is is 
repulsive for all values of 4 > and for all values of A < 2. 

In conclusion, we have found explicit formulae for the /c-string tension in both the weak 
and strong coupling phases of 2 dimensional lattice Yang-Mills theory. These formulae differ 
from the continuum result, which is the quadratic Casimir of the representation. They 
nevertheless share the property that the interaction between elementary strings is attractive 
when the quark sources are in the antisymmetric representation and repulsive when they are 
in the symmetric representation. Furthermore, in the weak coupling phase the k — > N — k 
duality is realized by a first order phase transition at k = y. An interesting check of the 
validity of these results is to examine the continuum limit. The weak coupling phase has a 
good continuum limit as A — > 0. In that limit, we should send the lattice spacing to zero, 
a — >■ while holding the string tension fixed. If we require that the k = 1 string tension in 
inverse distance units 

1 4 

0i = — In r (33) 



a 2 4- A 

remains finite, we must tune the coupling constant as A 
/c-string tension fl3"Ul) . we obtain 



4 - 4e~ 



Using this in the 



Similarly, for the symmetric representation in ( 132]) . the continuum limit yields 
a Sk = k (1 + A) ai - |L ((1 + *) aVi)Vi + 0(a e ) . 



(34) 



(35) 



The leading terms are in agreement with the Casimir scaling ()3]) which is known to be the 
exact behavior of continuum 2-dimensional Yang-Mills theory. 
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